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Abstract 

We derive a general WKB energy splitting formula in a double-well potential 
by incorporating both phase loss and anharmonicity effect in the usual WKB 
approximation. A bare application of the phase loss approach to the usual 
WKB method gives better results only for large separation between two po¬ 
tential minima. In the range of substantial tunneling, however, the phase loss 
approach with anharmonicity effect considered leads to a great improvement 

on the accuracy of the WKB approximation. 

03.65.Sq 


Typeset using REVTgX 


1 


I. INTRODUCTION 


The quantum tunneling in the double-well potential, V{x) = \x^ — kx'^{k,X > 0), is a 
longstanding and well-known problem in quantum mechanics. One of the main interests in 
this problem is how accurately one can calculate the energy splitting of degenerate energy 
levels due to quantum tunneling. The nonperturbative nature of the double-well potential, 
which stems from the quartic term, does not allow standard perturbation treatment. Instead 
several alternative methods have been proposed. Among those there are the instanton 
method , the WKB approximation , and a numerical method ^ . 

In the instanton method, one uses imaginary time path integral to obtain the classical 
action - the so-called euclidean action. Qualitatively, this method is useful in understanding 
the quantum tunneling which has no classical correspondence. However, the quantitative 
calculation of the energy splitting in the double-well potential by this method is inaccurate 
because the euclidean propagator can be obtained only in the limit of infinite separation 
between the two potential minima [Q, which corresponds to zero tunneling probability. 
Thus, the validity of the instanton approach is restricted to the case of which the two 
potential minima are far apart, and its accuracy is expected to be reduced as they become 
close. 

In the case of WKB method, one uses the semiclassical approximation in which the de 
Broglie wavelength \{x) = h/— V (x)) is assumed to be very short and varies suffi¬ 
ciently slowly. The main defect of this approximation comes from the inevitable difficulties 
in treating the divergence of A(x) at a classical turning point. Recently, a phase loss ap¬ 
proach has been proposed to improve this shortcoming for general barrier potentials . In 
the usual WKB approximation, based on the short wavelength condition, the phase change 
of the WKB wave function at a classical turning point is taken to be 7r/2. However, since the 
wavelength should not be arbitrarily short, the phase change should become a nonintegral 
multiple of 7r/2. Based on this idea the authors of Ref.5 have shown that the accuracy of 
the WKB wave function is greatly improved. 
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In this paper we use this phase loss approach to calculate the WKB energy splitting 
in the double-well potential. Although the authors of Ref.5 applied this approach to sev¬ 
eral barrier potentials, the double-well potential, which is one of the prototypes in tunneling 
problem, was not treated. In addition to this we also consider the anharmonicity effect of the 
double-well potential. In most previous works @,|,§| the energy eigenvalue of a unperturbed 
harmonic oscillator has been used for the calculation of the energy splitting, which corre¬ 
sponds to the case of the infinite separation of the two potential minima, and hence exposes 
same defect as the instanton method. The present work, however, incorporates anharmonic¬ 
ity into the WKB formalism, which gives a more realistic model, and consequently a more 
improved result in the region over which the tunneling is appreciable. By using the phase 
loss approach with the anharmonicity effect we will present a general WKB energy splitting 
formula in a double-well potential. From this, the usual WKB energy splitting formula can 
be obtained by both taking short wavelength limit and using energy level corresponding 
to the inhnite separation of the potential minima. Also, if we take the short wavelength 
limit only we can derive an anharmonic WKB energy splitting formula. We then present 
comparisons between numerical results obtained from these formula ; for completeness, we 
also include results from the instanton method. We find that a bare application of the phase 
loss approach to the usual WKB method leads to better results only in the limits of large 
separation. If we consider the anharmonicity effect, however, the phase loss approach greatly 
improves the accuracy of the usual WKB approximation, especially in the region where the 
tunneling is appreciable. In the following section, we will first derive a general WKB energy 
splitting formula. In section III, we will show comparisons between results of various cases 
and give discussions. Finally, there will be a summary in section IV. 

II. ENERGY SPLITTING FORMULA 

We consider a double-well potential of the type given in the introduction, which is shown 
in Fig.l. The Hamiltonian for this potential is 
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H = + \x‘^ — kx^ 


( 1 ) 


in which p is the particle momentum, and we have taken the particle mass m = 1/2. We will 
also take ^ = 1 throughout this paper. In the WKB approximation we have phase integrals 
in the two classically allowed regions: 

n r-i 

a= Pn{x)dx, b= Pn{x)dx, (2) 

J a J —o. 

and inside the central barrier we also have 

/ a 

Kn{x)dx, (3) 

-a. 

where 

Pn{x) = ^jEn-V {x) = iKn{x) (4) 

with En being an energy level in one of the wells. When we include the anharmonicity effect 
this energy level becomes, to second order correction, 

En = (5) 


where 


Bf = (2n + IjVa - 


( 6 ) 


is the energy level of a unperturbed harmonic oscillator and is the second order correction 
term given as 


. 3A, 2 U A2 

dn = —+ n -\ —)- 

k ^ 3 ^ QAk^^M 


(34n^ + 51n2 + 59n + 21). 


( 7 ) 


To obtain the WKB energy splitting formula in the double-well potential we follow the same 
method as described in Ref.7 except using arbitrary phase losses (/o(= 0_a), 0.y(= ip-'y) at 
each turning point instead of 7r/2 which corresponds to the short wavelength limit. After 
long but straightforward calculations we obtain 
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1 4e‘^tan^ + 0^) + 4e'^ 
T tan^ \ {.4>a + ^ 7 ) — 1 


( 8 ) 


where 


T = 


(9a 


1 n dx 


(9) 


^(9i9„/ 2 Ja Pn{x) 

To find the phase losses at the turning points we modify the matching condition suggested 
in Ref.5. Since each well is asymmetric the phase changes at a and 7 are different. At 
X = ±7 the phase losses are given by 

-1 r 1 ^'( 7 ) 


= —2 tan 


Pmln V’(7) 


( 10 ) 


where 


Pmln = \jE„- V{±tJk/2\), 


( 11 ) 


and 



(12) 


(13) 


The phase losses <Pa) at ± 0 ; can be obtained from the following algebraic equation : 


Pmin tSiIl I Q 


(14) 


^(-«) .. 2 ^ 

In this case the Lippmann-Schwinger equations should be written by 

/a p'y 

{x + a)K^J‘^{x)^i{x)dx + {x + a)p^^‘^{x)'ilj 2 {x)dx 

-a Ja 

roo 

+ / {x + a)Kl[‘^{x)%lj 2 ,{x)dx (15) 

J 7 

'ilj\-a)= f Kl^‘^{x)^jJl{x)dx + ( p^J‘^'ilj2{x)dx + f K^J‘^{x)'ilj‘i{x)dx, (16) 

J — Oi J Oi J'y 

where ipiix), iIj 2 {x), ipsix) are the WKB wavefunctions inside the central barrier, in the right 
well, and in the right-barrier, respectively. By applying the WKB connection formula at 
each turning point we can obtain 
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'iplix) 

'lp2{x) 

^3{x) 


= e 


cos A - 2 sin A 

““^cosA sin — J pn(x')dx'^ — 46 “^ sin A cos 


Pn{x')dx' 


-e“^ + 
,4 


sm 2 A 


( 17 ) 

(18) 

(19) 


with 


A = a — 


"h (f)'y 


We note that Eq.(8) is a general WKB energy splitting formula with both anharmonicity 
effect and phase loss. In this formula, if we use the unperturbed energy level we have 
the WKB formula with phase loss only: 

1 4e'^^°'tan^ |(0Q, + 0..,,) + 4e'^^°' 


where 


with 


A z? — 

^^-p-T(o) 4e25(°)tan2 1(0, + 0.,)-l ’ 


y(o) _ 1 r dx 
2 


/ (o1^ K ^ I ^n\x)dx 

Pn [X) 


( 20 ) 


pi"> = i/rf’ - V(x). 


Also, If we let 0^ = 0.y = 7r/2 in Eq.(8), it reduces to 

AEan = (21) 

which corresponds to the case of the WKB approximation with anharmonicity effect only. 
Finally, by both taking 0 q, = 0.^ = 7 r/2 and using Eq.(20) the usual WKB energy splitting 
formula can be obtained to be 




^ s(o) 
T(o) 


( 22 ) 
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III. NUMERICAL RESULTS AND COMPARISONS 


To examine whether and how much the general WKB energy splitting formula improves 
the WKB approximation we calculate the various energy splittings derived in the previous 
section numerically and compare each numerical result with the exact values obtained by 
the nonperturbative method [4], For completeness, we also include numerical results from 
the instanton method by which the energy splitting is obtained to be [2] 


A^in 


ttA 


kVn 

e 3a . 


(23) 


The numerical calculations have been performed for n = 0 (the ground state) case with 
k = 1, and we hnd AE^s for different A’s ; the n = 1 case does not give enough values for 
comparisons. For direct comparisons, we compute the ratios of each numerical result and 
the exact value taken from Ref.4, then hnd the respective differences of the ratios from 1. 
The hnal results are listed in Table 1. 

An immediate observation from this table is, as remarked in the introduction, that both 
the instanton (the hfth column) and the usual WKB (the hrst column) methods give good 
estimate only for small values of A which correspond to the limiting case of inhnite separation 
of the two potential minima, and their results are getting worse as A increases IQ]. It is also 
clear that the usual WKB approximation is better than the instanton method, which is 
consistent with the result of previous work [4], 

For the presence of tunneling in the double-well potential should the separation between 
the two potential minima be neither too small nor inhnitely large. When the two minima 
are too close the height of the barrier becomes lower than the lowest energy level in each 
well. This can be seen from the first four columns where energy splittings for A = 0.2 do 
not appear. In the opposite limit of A —0 the two potential wells become two independent 
harmonic oscillator potentials, which does not allow tunneling. In order for the tunneling 
probability not to vanish the two potential wells should not be inhnitely apart, which requires 
the consideration of the anharmonicity effect in calculating the energy splittings. From 
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Eq.(5) and Eq.(7) this effect is seen to be little for small values of A. The calculation of the 
energy splitting due to tunneling in this region is, then, governed by the usual WKB formula. 
As A increases, however, the two wells approach each other, and the anharmonicity effect 
become more important. Especially, for an appreciable tunneling probabiliy the anharmonic 
property of the double-well potential should be incorporated in the calculation of the energy 
splitting. Thus, the result from the WKB formula with anharmonicity effect is expected 
to be improved as A increases. We can see this behavior from the comparisons of the first 
column and the third one (the WKB result with anharmonicity only). For A < 0.05 the 
usual WKB results are more accurate than that with anharmonicity effect. But, for A > 0.07 
the results from the WKB formula with anharmonicity effect are much more accurate than 
the usual WKB result. Note that for A = 0.15 the error of the anharmonic result is less than 
2%, while that of the usual WKB result has almost 20%. This illuminates that the usual 
WKB approximation becomes unreliable for small separation distance. 

We now compare the first column with the second one. For small values of A(< 0.03), 
the phase loss approach is better than the usual WKB estimate. For A > 0.04, however, the 
result from the phase loss is more inaccurate than the usual WKB result. By comparing 
the third column with the forth one we also find that the results from the general energy 
splitting formula which includes both the phase loss and the anharmonicity effect is slightly 
better than those from the WKB method with anharmonicity effect only. According to Ref.5 
the calculation of phase loss at turning point is based on the assumption that the WKB 
wavefunction inside a well is proportional to the exact wavefunction which can be obtained 
from Schrodinger equation around the potential minimum. The condition for the validity of 
this approximation is that the dimensionless parameter Pmina^o, where xq being the distance 
between the point of potential minimum and the turning point, should be small. In the 
case of usual WKB method, while Xq increases with A, Pmin = \[e^ is independent of A 
so that PminXo becomes larger with increasing A. Thus, the application of the phase loss 
approach to this case is expected to be valid only for small \'s. On the other hand, when 
the anharmonicity effect is included, from Eq.(ll) decreases with A. Although the Xq 



value in anharmonic case also increases with A the dimensionless parameter Pmina^o remains 
small even for large values of A. Consequently, the incorporation of the anharmonicity effect 
in the WKB approximation makes the phase loss approach valid even for large values of A. 
Therefore, the entries in the second and the forth columns are reasonable results. 

Finally, we discuss the occurence of a larger energy splitting in phase loss approach than 
in the usual (or anharmonic) WKB approximation. If we compare the hrst column (or 
the third column) with the second one (or the forth one) we can notice that the energy 
splitting Ai^anp (or AE^p) obtained by the phase loss approach is larger than AF'an (or 
AE^j). The reason for this can be explained as follows. Since the usual WKB approximation 
is based on the condition of short wavelength limit its formalsim is more particlelike. Thus, 
in this semiclassical limit, the tunneling effect, which originates from pure wave property, 
is ill treated. However, in the phase loss approach we use the WKB method away from the 
short wavelength limit, and hence the approach is more wavelike. Therefore, the phase loss 
approach treats the tunneling effect better, by which a larger energy splitting follows. 

IV. SUMMARY 

We have investigated the phase loss approach in the WKB energy splitting formula about 
a double-well potential. By incorporating the phase loss approach with the anharmonicity 
effect in the usual WKB approximation we have obtained a general WKB energy splitting 
formula. A bare application of the phase loss to the usual WKB approximation yields a 
slight improvement only in the limits of large separation of the two potential minima. When 
the anharmonicity effect is included, however, the phase loss approach signihcantly improves 
the accuracy of the usual WKB method, especially in the range of appreciable tunneling. 

In applying the phase loss approach to the WKB approximation it is crucial to calcu¬ 
late the phase change at turning point accurately, for which there exists a condition such 
that the dimensionless parameter PmmXo is small. As noted in Sec.Ill the inclusion of the 
anharmonicity keeps this condition over a region where otherwise the condition becomes 
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invalid, so that a correct phase change at turning point can be obtained. It is thus impor¬ 
tant to consider the anharmonicity effect in the application of the phase loss approach to an 
anharmonic potential such as double-well potential. 

Finally, We have also found that the magnitude of energy splitting obtained from phase 
loss approach is larger than that of the WKB energy splitting obtained in the limit of short 
wavelength. 
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FIGURES 

FIG. 1. One-dimensional double-well potential. E is the ground state enrgy, and ±a, ±7 are 


the classical turning points corresponding to E. 
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TABLES 

TABLE I. Ratios of numerical result and the exact value taken from Ref.4. All entries are the 


differences of the ratios from 1. 


A 

A^’u 2 

A^/exact 

A£^up 

A£/exact 

A£^an 

A£/exact 

A£^anp 2 

A^/exact 

A£/exact 

0.02 

-0.03075 

-0.03046 

-0.07359 

-0.07357 

0.04524 

0.03 

-0.01297 

-0.01227 

-0.07382 

-0.07376 

0.07020 

0.04 

0.00543 

0.00676 

-0.07263 

-0.07252 

0.09773 

0.05 

0.02466 

0.02691 

-0.07006 

-0.06990 

0.12809 

0.07 

0.07348 

0.07883 

-0.05426 

-0.05396 

0.20694 

0.1 

0.15125 

0.16762 

-0.02278 

-0.02157 

0.35392 

0.15 

0.21684 

0.25018 

0.01799 

0.01829 

0.66831 

0.17 

0.07471 

0.08135 

-0.00882 

-0.00804 

0.77935 

0.2 

none 

none 

none 

none 

0.91421 
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